STABILITY OF THE GAUGE EQUIVALENT CLASSES IN 
INVERSE STATIONARY TRANSPORT IN REFRACTIVE MEDIA 



STEPHEN MCDOWALL, PLAMEN STEFANOV, AND ALEXANDRU TAMASAN 

Abstract. In the inverse stationary transport problem through anisotropic atten- 
uating, scattering, and refractive media, the albedo operator stably determines the 
gauge equivalent class of the attenuation and scattering coefficients. 



1. Introduction 

This paper concerns the problem of recovering the absorption and scattering prop- 
erties of a refractive medium from boundary knowledge of the albedo operator. The 
medium M C M. n , n > 2, is a bounded domain with smooth boundary, endowed with 
a known Riemannian metric g. The free moving particles travel through M along 
the geodesies. In the stationary case the propagation of particles is modeled by the 
linear transport equation 

(1) — Vu(x, v) — a(x, v )u(x, v) + J k(x, v', v)u(x, v ') du x (v ') = 0. 

J S X M 

In the equation above u(x, v ) denotes the density of particles at position x with 
velocity v in S X M, the unit tangent sphere at x. The operator T> is the derivative 
along the geodesic flow: For a given point (x,v) G S X M, if j( x ,v)(') denotes the 
geodesic starting at 7( Xi „)(0) = x with initial velocity jr x>v )(0) = v, then 

d 

(2) Vu(x,v):=— u(7(x,«)(<)), 

where, for brevity, we use the notation j( X}V )(t) = (j(x,v)(t),j(x,v)(t))- If g is Euclidean 
then T> is the directional derivative: T>u(x,v) = v ■ V x u(x,v). The measure du x (v') 
in ([1]) is the volume form on S X M induced from the volume form on T X M (the 
tangent space to M at x) determined by g at x. The resulting (Liouville) form on 
SM is preserved under the geodesic flow of g, see [22]. The attenuation coefficient 
a(x,v) in ([T]) quantifies the rate at which particles are lost from the point (x,v) in 
phase space due to absorption and scattering into new directions. The scattering 
coefficient k(x,v',v) represents the probability that a particle at position x with 
velocity v' G S X M will scatter to have new velocity v G S X M. 

The boundary measurements are described by the albedo operator A: Let F± = 
d±SM = {(x,v) G dSM : ±{v,v x ) > 0} denote the "incoming" and "outgoing" 
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bundles, where v x is the unit outer normal vector to the boundary dM at x and (■, ■) 
is the inner product, each with respect to g at x. The medium is probed with the 
given radiation 

(3) u\ v _ = it_, 

and the exiting radiation is detected on r + . The albedo operator takes the incoming 
flux to the outgoing flux at the boundary: Ail- = u|r+- 

The inverse boundary value problem is to determine the coefficients a and k from 
the knowledge of A. 

When the attenuation a is isotropic (v- independent), there is a large collection 
of uniqueness results under varying assumptions on the parameters; see j2] for a 
comprehensive account, and [20] for a recent survey of numerical methods. The 
works below are based on the singular decomposition of the Schwartz kernel of A, an 
idea first introduced in [5] and [7J; see also [5]. In the Euclidean setting, uniqueness 
of a(x) (dimensions two and above) and k(x, v', v) (dimensions three and above) was 
proven in [7J under some minimal restrictions guaranteeing only that the forward 
problem is well-posed. For sufficiently small k this result was extended to dimension 
two in [IS]; see also [H [27]. Stability results are proven in [Ml HOl HH [28] with 
the most general result (in Euclidean geometry) in [3]. When no angular resolution 
is measured in the outgoing flux, the singular decomposition of the new boundary 
operator has been used to recover an isotropic coefficient a and the spatial part of k 

m [nam. 

The case of a Euclidean metric corresponds to transport in materials with a con- 
stant index of refraction. If the index of refraction is isotropic, but varying, then 
(PQ) can be derived as a limiting case of Maxwell's equations with non-constant (but 
isotropic) permeability, resulting in a metric which is conformal to the Euclidean 
metric ([I]). For a general metric, we consider |1| as a model for transport in a 
medium with varying, anisotropic index of refraction. When the attenuation is as- 
sumed isotropic, uniqueness results in Euclidean geometry are extended to the Rie- 
mannian metric setting in [Til C31 O 03] ■ There and here, the manifold is assumed 
to be simple as follows. 

Definition: (M, g) is called simple if it is strictly convex, and for any x G M the 
exponential map exp x : exp~ x (M) — > M is a diffeomorphism. If M is two dimen- 
sional we have the following additional assumption: Let k be the maximum sectional 
curvature of M. If k> , then we also assume diam(M) < 7r/i/«. 

The works mentioned above concern the media with an isotropic attenuation char- 
acter. However, since the attenuation is a combination of absorption and loss of 
particles due to scattering: a(x, v) = a(x, v) + J s M k(x, v, v') du x (v'), even when the 
absorption part is isotropic (er = a(x)), if k depends on two independent directions 
the resulting attenuation is anisotropic. Evidence of anisotropy in biological tissue 
has been observed experimentally, see |12j . 

When the attenuation coefficient is anisotropic, it is possible to have media of dif- 
fering attenuation and scattering properties which yield the same albedo operator. 
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Moreover the non- uniqueness is characterized by the action of a gauge transforma- 
tion [23]: see (jl]) below. The same algebraic structure of non-uniqueness is valid in 
refractive media [T6j . 

In Theorems 14.11 and 14.21 we show the stability of the gauge equivalent classes 
occurring in refractive media, thus extending the results from the Euclidean case 
in [T7]. We also generalize a stability result in [3] from Euclidean to Riemannian 
geometry. 

The algebraic structure of non-uniqueness in [23l 116] can be readily observed. 
Indeed, if G L°°(SM) is positive with 1/0 G L°°(SM), V<p G L°°(SM) and such 
that = 1 on dSM. Set 

/,\ -/ \ / \ ,, \ 7/ / \ k(x,v',v)(j)(x,v) 

(4) a{x,v) — a{x,v) — V log (p(x,v), k{x, v , v) — — 



[x, v') 



Then u satisfies (OQ) if and only if u = 4>u solves 

— Vu(x, v ) — d(x, v )u(x, v ) + / k(x, v ', v )u(x, v') du x (v ') = 0. 

J S X M 

Since = 1 on T, u = £t there, so the albedo operator ^4 for the parameters (a, fc) is 
indistinguishable from the albedo operator A for the pair (a,k), i.e. A = A. This 
motivates the following definition in [23| [T6] . 



Definition 1.1. Two pairs o/ coefficients (a, fc) and (a, fc) are called gauge equivalent 
if there exists a positive map G L°°(SM) with 1/0 G L°°(SM), V<j> G L°°(SM), 
and = 1 on T, snc/i inat (j4j) holds. We denote this equivalence by (a, k) ~ (a, fc). 

The relation defined above is reflexive since (a, k) ~ (a, via = 1; it is symmetric 
since (a,k) ~ (d,k) via yields (a, A;) ~ (a,k) via 1/0; and it is transitive since 
if (a, k) ~ (a, fc) via and (a, fc) ~ (a', fc') via then (a, fc) ~ (a', k') via 00. 
Therefore one has the multiplicative group of gauges acting transitively (since any 
equivalent pair are related by some gauge [231 [IB] ) on the equivalent class of a pair 
of coefficients. We denote the equivalence class of (a, k) by (a, k). 

2. Transport of the data to a larger domain 

Due to the method of proof, the total travel time of each particle in M has to be 
uniformly bounded away from zero. This can be done without loss of generality by 
doing the measurements away from the boundary DM. More precisely, let Mq be a 
slightly larger domain strictly containing M. The metric g can be extended to go on 
M in such a way that (M ,g ) still remains simple [25] . 

As in [T7], we reduce the problem in M to one in Mq: Let (a,k) and (a,k) be 
coefficients for which the forward problems in (M, g) are well-posed, and A and A 
be their corresponding albedo operators. Defining a = d = k = k = 0in Mq \ M, the 
forward problems in (Mo, go) are also well-posed and the albedo operators Aq and 
.Ao are we ll defined maps between functions on 

:= d±SM = {(x, v) G 8SM : ±{v, v x ) > 0} 
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(now v x is the outer unit normal vector to dM at x, with respect to the extended 
metric go). As in [T7], when the two forward problems for M are well-posed in L p , 
1 < p < oo, the following isometric property holds: 

(5) ||*4 - A\\c(LP(r^,dn);LP(r + ,dn)) = HA) - AoWc^^^^LP^dno))- 

The measure d\i (and, analogously, dfio) in (jSJ) is defined as follows: Let dY? n ~ 2 be 
the volume form on Y± obtained by the natural restriction of the volume form on 
SM to r±. Then, by extending dH 2n ~ 2 as a homogeneous form of order n — 1 in \v'\, 
we have that (iE 2n_2 (x', v') coincides with the volume form on SM; see [22] for 
details. We define 

(6) dft(x',v') = \(v',u x ,)\dJ: n ~ 2 (x',v'). 

Remark: The proof of (j3J) is essentially identical to that in the Euclidean case 
as presented in [T7] due to the following invariant property of the form dfi: Fix 
(x' ,v' ) G r ± . Let dCl be any surface so that the geodesic issued from (x' ,v' ) hits it 
transversally. Then the geodesic flow defines a natural local "projection" near (x' , v' Q ), 
of r± onto T±. Let djx be the measure on r± defined analogously to (jSJ). Then the 
pull back of dfl is d\x (see the proof of Lemma 4.2.2 in [22]). Given (x,v) G SM, we 
define the "forward/backward travel time to the boundary" functions 

t±(x,v) = mm{t > : ^ x>v) (±t) G <9M }, 

and let t(x,v) = t + (x,v) + T-(x,v) be the total travel time of the free particle (of 
velocity v at x) through M . Since dist 50 (M, dM ) > 0, we have 

(7) c := inf {r{x, v) : (x, v) G ~SM} > 0. 

Using the isometry property we can considered the inverse problem in the larger 
domain (Mo, go) with the albedo operators now acting between T±. Equivalently, for 
the original problem in (M, g) we may work without loss of generality with coefficients 
a, k of (a priori fixed) compact support in M. 

To simplify notation, while still working in the larger domain, we drop the index 
throughout the remaining of the paper: thus M becomes M, T± becomes r±, etc. 

3. The singular structure of the albedo operator's kernel 

In this section we recall the singular decomposition of the Schwartz kernel of the 
albedo operator for the two cases separated by dimension. 
We work within the class of admissible coefficients: For n > 3 

(8) (a, k) G L°°(SM) x L°°{SM, L\S X M)), 
and for n = 2 

(9) (a, k) G L°°{SM) x L°°(S 2 M), 

where S 2 M := {(x, v', v) : x G M, v', v G S X M}. Note that the gauge transformations 
(J4]) preserve the admissible classes in (jSJ) and ([9]). 



Stability of the gauge equivalent classes 



5 



Moreover, either one of the following subcritical conditions that yield well-posedness 
for the boundary value problem (TjQ) and Q is assumed to hold: 

(10) ess sup (XjU)e5M 
or 

(11) a(x,v) — / k(x,v,v f ) duj x {v') > 0, a.e. (x,v) G SM; 

Js x m 

see, e.g., [3 EH £9] . 

The right hand side of ([T]) defines a closed, unbounded operator on L l (SM R ) with 
the domain {u G L\SM R ) : Vu G L l {SM R ), u\ r _ = 0}; see [HQS]. 

Proposition 13.11 below, which describes the terms in the expansion of the kernel 
of A, is proven in [13], the Euclidean equivalent appearing in [7]. We denote by 
5{x',v'}{%, v ) the delta-distribution on T + with respect to the measure d[i defined by 

/ ip(x, v)8{ x >yy(x : v) d/i(x, v) = tp(x', v') : Lp e c^°(r + ). 

Similarly, 5{ x }(y) is the delta distribution on M supported at x. 

If x, y e M let v(x, y) G S X M denote the tangent vector at x of the unit speed 
geodesic joining x to y (uniquely defined since M is simple), and d(x,y) be the Rie- 
mannian distance between x and y. Denote the total attenuation along the geodesic 
from x to y by 

(12) E(x,y) := exp|-y a(%^ X) y))(t)) dtj. 

Note that j( y ,v(y,x)) (d(y, x) — s) = —j( x ,v(x,y))(s), so when a depends on direction, 
E(x,y) ^ E(y,x). 

Proposition 3.1. |13] Let (M,g) be a smooth simple Riemannian manifold of di- 
mension n > 3. Assume that (a, k) are admissible and subcritical so that the forward 
problem is well-posed. Then the albedo operator A : L 1 (r„,<i/i) — > L 1 (T + ,dfi) is 
bounded and its Schwartz kernel a(x,v,x',v'), considered as a distribution on T + 
parameterized by (x',v') G T_, has the expansion a = a + ot\ + a 2 , where 

(13) a = E(j {x>v) (-T-(x,v)),x)5tf (xfy) ( T{x , >vl))} (x,v), 

/•T+(x',v') PT—(x,v) 

(14) ai = / E(y(s),x)E{x\z{t))k{z(t) 1 z{t),y(s))6 {y{s)} {z{t))dsdt 

Jo Jo 

y(s) = l(x,v)(s - T-(x,v)), z(t) = 7( x >y)(t), 

(15) a 2 GL 0O (r_;L 1 (r +J d//)). 

We note that k(z(t), z(t),y(s)) is only defined on the support of the integrand, 
namely when y(s) = z(r). 

When n = 2 the left hand side of ([1]) defines an unbounded operator of domain 
{u G L°°(SM) : Vu G L°°(SM), u\ r _ = 0}. Provided that © holds, and we have 



t(x,v) / k{x, v , v ') du x (v ') 



S X M 



< 1, 
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subcriticality f fTUj) or (fTTT) . this operator has a bounded inverse in L°°(S'M), see [I] 
and the singular decomposition of the albedo kernel is more explicit as follows. 

Given (x,v,x',v') G T + x T_, define x '■ T+ x T_ — >• {0,1} by x{ x i v ■> x ' ■. v> ) = 1 
if there exist < s = s(x,v,x',v') < T-(x,v) and < t = t(x,v,x',v') < t + (x',v') 
such that 7(a;,u)(s — r_(x,t>)) = j^y^t) (i.e., the geodesies intersect in M), and 
x(x, v, x', v') = otherwise. When = 1, let ip(x, v, x', v') be the angle 

between the tangent vectors of these geodesies at the point of intersection. 

Proposition 3.2. |l3j Let (M,g) be a two dimensional simple Riemannian mani- 
fold. Assume that (a, k) are admissible and that (Q holds. Then the albedo operator 
A : L°°(r_, dfi) — > L°°(r + , dfi) is bounded and its Schwartz kernel a(x, v, x', v'), con- 
sidered as a distribution on Y + parameterized by (x',v') G T_, has the expansion 
a = «o + ot\ + «2, where 

a = E(7( X7v) (-T-{x,v)),x)5tf (x , y) ( T ( x/>v , m (x,v), 

k(%',v')(t),%,v)(s - T-{X,V))) 



sm(ip(x, v, x', v'))\ 



a>i = x(x,v,x', v')E(x' ^( x >y)(t))E(-f^y)(t), x)J - 

0<a 2X < C\\k\\U {s2M) {l + log |sin( ^ (x ^ /;t/))| ). 

Here, J = J{x, v, x', v') is a function uniformly bounded < mi < J < m 2 < oo on 
T+ x T_ (see [HJ Proposition A}). 

4. Statement of the main results 

Let (B a , || • || Bq ) and (B k , \\ ■ \\ Bk ) be Banach spaces in which the attenuation and, 
respectively, the scattering kernel are considered, (a, k), (a, k) G B a xB k . The distance 
A between equivalence classes with respect to B a x B k is defined by the infimum of 
the distances between all possible pairs of representatives: 

A((a,k),(a,k)) : = inf max{||a' - a'|| Bo , \\k' - k'\\ Bk }- 

(a',k')€{a,k), (a' ,k')e{a,k) 

The following norms are used throughout 

Halloo = ess s\xp (x>v)eSM \a(x,v)\ : 

\\k\\oo,i = ess swp iX)V , )eSM / \k(x,v',v)\ du x (v), 

JS X M 

s ^P(x,v',v)esm\ k ( x , v ', v )\, 



, = ess 

ll^lli = / / / \k(x, v', v )| dx du x (v') du x (v ). 
Jm J s x m Js x m 

Case n > 3: Define the class 

Ut,, p '■= {{a, k) as in (JED : ||a||oo < E, ||/c||oo,i < p}- 

Theorem 4.1. Let (M,g) be simple Riemannian manifold of dimension n>3. Let 
(a, k), (a, k) G U^, p be such that the corresponding forward problems are well posed. 
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Then 

A((a, k), (a, k)) < C\\A - A\\c(L^(r-,dfiy,L 1 (r + ,d li )), 

where A is with respect to L°°(SM) x L l (S 2 M), and C is a constant depending only 
on S, p, n, and c in (J7]). More precisely, there exists a representative (a', /c') G (a, A;) 
swc/i that 

(16) ||a' — a||oo < CM ~~ ^■IU(L 1 (r_,dAt);L 1 (r + ,d/i)), 

( 17 ) W ~ Hi < CM - ^|U(Li(r_,dA t );i 1 (r+,^))- 

Case n = 2: From Proposition 13.21 above recall the Schwartz kernel of the albedo 
operator in the form 

a = A (x, v)5tf {xl v , ){T{x iy ))} (x, v)+/3 

where 

A (x,v) = E( 1{x>v) {-t-{x,v)),x) e L°°(T + ) 
P(x t v 7 rf,v')x\&nil>(x,v,a/,v')\ G L°°(T^ x T^). 

We define 

(18) \\A\U = max{|LA |U \\/3 X \ sin^Hloc}- 

By using the Remark in Section [21 the proof in [17] carries through verbatim to show 
that H.4. — *4.||* is preserved when transported from the boundary of the inner domain 
M to the boundary of the larger domain. 
Define the class 

Vt,, p '= {(a, k) as in @JJ) : Halloo < S, \\k\loo < p). 

Theorem 4.2. Let (M,g) be a two dimensional simple Riemannian manifold. For 
any E > there exists p > depending only on £ and (M, g) such that the following 
holds: if(a,k), (a, k) G Vs tP then 

A({a,k),{a,k)) < C\\A- A\U 

where A is with respect to L°°(SM) x L°°(S 2 M) and C is a constant depending only 
on E and (M, g). 

Note that p sufficiently small already yields a subcritical regime as in (TTOl) . 

5. Preliminary estimates 

In this section we extend a result from the Euclidean to Riemannian setting; see 
[31 Theorem 3.2] for contrast. In addition, the proof below allows for discontinuous 
coefficients, which is needed when transporting the albedo operator to the larger 
domain (and no boundary knowledge of the coefficients is available). 
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Lemma 5.1. There is a family of maps (/> E)X i v > G L 1 (T_,dfi), for (x' ,v' ) G T_ and 
e > 0, such that ||0 e ,x' o ,,,Ji,i(r_,<2/i) = 1 and, for any f G L°°(r_, djj) given, 



(19) lim / (p e y o y o (x',v')f(x',v')dfi(x',v') = f(x' ,v' ) : 



whenever (x' ,v' ) is in the Lebesgue set of f . In particular, ( TT9~j) holds for almost 
every (x' , v' ) G r_. 

For a measurable function / on IR n , The Lebesgue set of / is 

L, := {x : lim^-^i / \f(y) - f(x)\ dy = 0} 

r^O \B r (x)\ J Br{x) 

where B r (x) is the ball of radius r centered at x, and where | ■ | denotes Lebesgue 
measure. We point out that if G L} oc (R n ), then \R n \L f \=0 (P Theorem 3.20]). 

Proof. For (x' , v' ) G T_ and £ > sufficiently small, let (x', v') : U x W C M™ -1 x 
R n_1 ->■ dSM with x'(C/) C dM be a coordinate chart near {x' ,v' Q ) = (x'(0), v'(0)). 
Let dY> 2n ~ 2 (x', v') = J~gZ du dw be the local coordinate expression for the volume 



element (see ([6])). For (x', v ') G T_, define 

1 1 

° ° \Wx',v') \ y/g_(x',v') 

where <p(u) = l/(u> n -i) for \u\ < 1, (/?(«) = for \u\ > 1, and <p e (w) = e~ n+1 <p(u/e). 
By we denoted the volume of the unit ball in R n ~ l . Then, for any e > 0, 

f ip £ (u)du = 1 and, by using fl6]), we obtain 

(Pe,x' ,v' V')f(x', V') d/J,(x', v') 

1 



^( U (x'))Ve(w(x', V'))f(x', V') ===== dE- 2 (x', t/) 



(p e (u)(p e (w) f (x' (u) , v'(u, w)) du dw. 

-2 

Apply the equality above to / = 1 to get \\(p e ,x' y \\L 1 (r + ,d^) = 1- The conclusion 
follows from the approximation of identity for LP maps, see, e.g., (9j Theorem 8.15]. 

□ 

Consider F : dM x SM -> R defined by 

(20) F(x', y, w) = £(x', 7(y,«o(T+(l/, «>))), 

with E as in (TT2|) . F(x',y,w) represents the total attenuation along the broken 
geodesies from x' ->■ y -> 7 (2/)IU )(t+(?/, w)). 

Let (a, fc), (a, fc) be admissible pairs as in dSJ). Recall that, after the extension of the 
domain, the coefficients have fixed compact support away from the boundary dM. 
All the operators bearing the tilde refer to (a, k) and are defined in a similar way to 
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the ones for (a, k), i.e., A is the albedo operator corresponding to (a, k). Recall that 
n is the dimension of the space. To simplify notation let 



\A-A\\ := \\A-A\\ 



C{L^{Y^,d^);L 1 (Y + ,d^))- 



Theorem 5.2. Let (a,k),(a,k) be as in ([8]). For almost every (x' ,v' ) G T_ the 
following estimates hold: For n > 2, 



(21) 



< \\A-A\\. 



For n>3, with y(t) = ^f( x ' ,v' )(t) , 

T +( x O' v 'o) 



\ k - k\(y(t),y(t),w)F(x' ,y(t),w) du y (w) dt 



Sy{t)M 



(22) 



< \\A-A\\ + \\F-F\ 



k(y(t),y(t),w) dujy(w) dt. 



Proof. Let (x' ,v' ) G T_ be arbitrarily fixed and let <p e ,x' ,v' G L l {TJ) be defined as 
in Corollary 15.11 To simplify the formulas, since (x' ,v' Q ) is fixed, in the following we 
drop this dependence from the notation <p e = (j> EjX i v ' . 

Let A = Aq + Ai + A2 be the decomposition of the albedo operator given by 



Aif(x,v) 



oti(x, v, x', v)f(x', v) dji{x' , v), i = 0,1, 2, 



where ccj, i = 0, 1, 2 are the Schwartz kernels in Proposition 13. 1[ 

Let (f> G L°°(r + ) with ||0||oo < 1- Since ||^ e ||ii(r_) = 1, the mapping properties of 
the albedo operator imply that 



(23) 



<J)(x, v)[A — A]4> e (x, v) dfi(x, v) 



< \\A-A\\. 



Next we evaluate each of the three terms in J T <f)(x,v)[A — A](f> £ (x,v) dfi(x,v) by 
using the decomposition in Proposition 13. II and Fubini's theorem. 
The first term is evaluated using the formula (fl~3f) : 

(f){x,v)[A - A }4> £ (X,V) dfi{x,v) = / 4>(^ {xly) (T + (x',v')))<f) £ {x',v') 



d[i(x', v'). 



Since the integrand above is in L°°(r„) by applying (1T91) . we get for almost every 



{x' ,v' )er_ 

I ((j))(x' ,v' ) := \iml ((f),e) 

E—>0 



(24) =<f>{l M) (T + (x' ,v' )))(e 



Jo 
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To evaluate the second term we use the formula (fTSj) and let y = y(x',v',t) 

1{x',v')(t)- 



h((j),e):= / <j)(x,v)[Ai - Ai](pe(x,v) d/i(x, v) 
'r+ 

, rT+(x',v') p 

4> £ (x',v') dfi{x',v')\ / (f)(% tW) (T+(y,w))) 



r_ 

x 



F(x', y, w)k(y, y, w) - F(x', y, w)k(y, y, w) 



dw dt > . 



Apply again ([19]) for the continuous integrand above to obtain for almost every 
(x' ,v' ) G T_: with y = y(x' ,v' ,t) 

h((f))(x' ,v' Q ) := lim h{<f>,e) 

e— >>0 

(j)(l {y , w )(T + (y,w))) 

(25) x [F(x' , y, w)k(y, y, w) - F(x' , y, w)k(y, y, w)] du y (w) dt, 

or / x (0) = 7i,i(0) + /i, 2 (0) with 

(26) 

PT + (X' ,V' ) P 

h,x((j>) = / / (j)(% tW )(T + (y,w)))F(xQ,y,w)(k-k)(y,y,w)duj y (w)dt, 

J0 JSyM 

(27) 

\hM\< / \F-F\(x' Q} y,w)k(y,y,w)duj y (w)dt. 

JO JSyM 

Consider the third term 

J 2 (0,e)= / ^(t,u)[^4 2 - .4 2 ]0 e (x,?;) d/i(x,v) 
r+ 

' 4> e (x', v') dfj,(x', V ')^J 4>{ x i v )( a 2 — 0.2) {x, v, x', v') dfj,(x, 

By f|T5l) . the map (x',v') (-)■ J" r+ 0(x, w)(a 2 — o: 2 )(:r, v, v') d/i(x,v) is in L°°(r_), 
and then, by (|T9|) . we get for almost every (x ,t> ) G r_ 

(28) I 2 {<f)){x' ,v' ) :=liml 2 {(f),e) = / 0(a;, v)(a 2 - a 2 )(x, v, x' Q , v' Q ) d/i(x, v). 

Jr+ 

The left hand side of ( 1231) has three terms. We move the third term to the right hand 
side (with absolute values) and take the limit with e — > to get 

(29) 1/0(0) + /! (0)104^) < WA-AW+hMMA), a.e. (x' ,v' ) G r_, 

for any G L°°(r+) with ||^>||oo = 1- 

We note that the negligible set on which the inequality above does not hold may 
depend on 0. We will consider a countable sequence of functions 0, and since the 
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countable union of negligible sets is negligible, the inequality ( 129]) holds almost every- 
where on T_, independently of the term in the sequence. This justifies the argument 
below for almost every (x' ,v' ) in T_. 

In ( )29|) . we shall choose two sequences of to conclude the two estimates of the 
lemma. First we show the estimate (l2Tj) by choosing m G L°°(r + ) which are 1 in 
a shrinking neighborhood of {xq,vq) := l(x' ,v' )( T +( x 'o, v 'o))- First, define (j) m (xo,v) to 
be the indicator function for the set {v G S Xo M : \\v — Vo\\ g ( xo ) < 1/m}, then extend 
0m by 

t x v \ _ |° if d d M{x,x ) > 1/m, 

\<p m (x,V{v;x,x )) if d dM {x,x ) < 1/m 

where V(v; x, xq) is the parallel transport of v G S X M from x to xo- Then ( )24|) gives 

loWm) — e ^o'V — e v '^o'V 

independently of m. From f )25|) we have lim^^co ii (0 m ) = since for any t, the 
support of 0(7( y ( t ) jW )(r + (2/,u>))) in G S y ^M shrinks to From (l2"gj) we also 

have lim m _ il00 I 2 (|0 m |) = 0, since the support shrinks to one point. We use here the 
corollary of (TT5]) that (a 2 — a 2 )(-, •, Xq, Wq) G L 1 (r + ), for a.e. (xq^o) G T_. 

Next we prove the estimate (1221 . Recall that now n > 3. For m > 0, let N^ x ' 0)V ' o ^ q C 
M be the tubular neighborhood of the geodesic y(t) = r Y( x , ,i/ )(t), < t < r + (x' ,v' ), 
of radius 1/m. We now define a sequence m G L°°(r + ). Set <fi m {x,v) = if 
z e N( x > otV > o ) 7 g; note that / o (0m) = for all m. For (x,u) G T + with x ^ N( x > o>v > o ) >q , 
4> q (x,v) = if the geodesic z(s) = 7( X]V )(s), — T_(x,t>) < s < does not intersect 
N( x ' 0>v ' ) >q . When 7(x,«)(-) does intersect N( x < o>v > o ) >q , let < t(x,w) < t + (x' ,v' q ) and 
— r_(x, v) < s(x,v) < be such that 

d g (y(t(x,v)),z(s(x,v))) = mm{d g (y(t),z(s)} 

and define 

(j) m (x,v) = sgn(&; - k)(y(t(x,v)),y(t(x,v)),V(z(s(x,v)); z(s(x,v)),y(t(x,v)))). 

Notice that when (x, v) is of the form 1(y(t),w){T+.{y{t), w)), w G S y ^M, that is (x' Q , v' Q ) 
and (x,v) are the beginning and end of a single-scattering broken geodesic, <fi m {x,v) 
takes the sign of k — k at the point of scattering. Note also that the support of m 
shrinks to a negligible set in T + as m — > oo since n > 3. 

Now apply the estimate fl29|) to m and use Jo (0m) = to get 

\h,i(4> m )\(x' ,v' ) < \\A- J[||+/ 2 (|0 m |)(x / o ,^) + |/ 1>2 (0 m )|(x / o ,^). 

Since the support of m shrinks to a set of measure zero in T + as m — > oo, we 
get for almost every (x' ,v' Q ) G T_, 

li m m— >oo 

/ 3 (|^ m |)(xo,^) = 0. Finally, noting 
that |/i,i(0 m )| = ii,i(0 m ) and applying (127)1 . from (|26|) we obtain for almost every 
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(4«o)er. 

r+{x' 0> v' ) 



\ k ~ M{y(t),y(t),w)E(x' ,y)E(y(t),'y(y( t)tW - ) (T + (y(t),w))) du y (w) dt 

-T + {x' ,v' ) 



lim I hl ((f> m ) < \\A-A\\ 



The estimate (1221) in the theorem follows. 



\ F - F \{x' ,y,w)k(y,y,w) duj y (w) dt. 



□ 



6. Stability modulo gauge transformations 

In this section we prove Theorem 14.11 

We start with two pairs (a, k), (a, k) G U^ p and let 

t := ll.A-.4ll. 



We shall find an intermediate pair (a', k') ~ (a, k) such that (fl6l) and ffTT|) hold. 
Define first the "trial" gauge transformation: 

(30) <p(x,v) := e -/o T " (a: "' ) (a-«)(7( a! ,v)( S -r-(^)))<i S) a e e £M_ 

Then 9? > 0, <^| r _ = 1, Vtp(x,v) G L°°{SM) and 

(31) 5(x, i>) = a(x, i>) — T>log<f(x, v). 

Note, however, that (p\r + is not equal to 1. We begin by estimating (fi\r + - By (j21 
we have for almost every (x' , v' ) G T_ 



< e. 



Changing variables t = t+(xq,v' ) — s and denoting (x ,Vq) = j( x i v n(r + (x' , v' )) we 
get 



(32) 



pT (3:0,1)0) . . .., r T — { x > v () ) / -> / / \\\ 

Jo a {l(x ,v )(t-T--{xo,v ))) dt _ - J a vy(a, ,v ) (bo,vo))) dt 



< £. 



When (x' Q ,v' Q ) covers T_ almost everywhere we get (xq,Vq) covers T + almost every- 
where. 

By the Mean Value theorem applied tou^ e~ u we obtain the lower bound 



(% ,v ){t~T-(xo,v ))) dt _ g-/ T iX °' V0) a(j( XOiVO )(t-T-(x ,v ))) dt 

T-(X ,V' ) 

(a - a)(% 0tVo )(t - t-(x ,vq))) dt = e~ u °\ \ogtp{x , v )\ 



(33) >e- diam W E |log^(xo,^o)| 

where u = u (x ,v ,a,d) is a value between the two integrals appearing in the 
exponents in the left hand side above, and <p is defined in fl30|) . 

From (132]) and (133]) we get the following estimate for the "trial" gauge ip: 

(34) I log <p(x, v) I < e diam{M)s e, a.e. (x, v) G T+. 
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The "trial" gauge <p is not good enough since it does not equal 1 on T + . We alter 
it to some (p £ L°°(SM) with T>\og<p £ L°°(SM) in such a way that (p\asM — 1- 
More precisely, for almost every (x,9) £ SM, we define <p(x,v) by 

(35) log<p(z,i;) := logy(x,t>) - - ' log<p(% iV) (T + (x,v))). 

Since < r_(x, v )/r(x, v ) < 1 we get (p £ L°°(SM), and clearly <p\osm = 1- Since 
Vt(x,v) = X>log^(7( a . jV )(r + (a;,'u))) = and Vt^(x,v) = 1, 

/o R N th -/ N t,, / \ lQg^(7(^)(r+(x,t;))) roo ^ ca ^ 

(36) X>log</?(x, = Vlog<p{x,v) — r -EL (SM). 

t(x, V) 

Define now the pair (a', k') in the equivalence class of (a, k) by 

(37) a'(x, v) := a(x, v) — T>log<p(x, v), 

(38) kf(x t J,v):=-?}^k{x t J,v). 

(pyx, v') 

Now A', the albedo operator corresponding to (a', k'), satisfies A' = A, and 

1 1 1^4. t^4. 1 1 — 1 1 ^4. 1 1 — • 

Next we compare the pairs (a', k') with (a, k) and show them to satisfy ([TBI) and 
(jTTJ) . Using the definitions (j3~T]) . ( l3"7|) . the relation f )36|) . and the estimate (}3"4"j) for </> 
on r + , we have for almost every (x, v) £ SM: 

\a{x, v) — a'(x, v) \ = | [5 — a](x, u) + [a — a'](x, v) | 

= \V log </?(x, v) — V log </?(x, v) | 

/o^ |iog v (tc^Ct+C^, «))) I , 

y6i>) — r S £ 

ryx,v) ryx,v) 

Since the coefficients are supported away from dM (by construction of M) such that 
(J?D holds, following f )39|) we obtain the estimate (ITB"]) in the form 

gdiam(M)S 

(40) ||a-a'||oo<e , 

c 

with Co from ((7j). 

Up to this point, all the arguments above also work for two dimensional domains. 
Next we prove the estimate (fT7|) . These arguments are specific to three or higher 
dimensions. Recall the formula ( 120]) adapted to a': let x' £ dM, y £ M and w £ S y M 
and let v' £ S x 'M,t > be such that y = j( x ',v')(t)- Then from (HOI) . 

gdiam(M)S gdiam(M)£ 

|a'(x, v)|<e — -, r — h E, and llo'lloo < £ h S, 

r(x, u) c 

so (using the first, and the fact that r is constant along geodesies), 

\F'(x',y,w)\ = e~J° a{ri ^'' v '^ s " e~J° a ^(y> w n s >> ds 

(41) > exp (_2( ee diam ( M ) s + diam(M)S) . 



e diam(M)E 
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Using the non- negativity of a and a 1 we estimate 
\[F-F'](x',y,w)\< 



e ~ So "■(%', v')( s )) ds _ g-/o lV ')(s)) ds 



< 



(42) 



[a - a'](%> )V >)(s)) ds 

t , r+(y,w) 
r(y, w) 



< £e diam(A/)E / * 

\t(x' v') 



r+{y,w) 

[a - a'}(% tW )(s)) ds 
< 2ee diam ^ M ^ s 



by ( 1391) . We now apply the lower bound for F' in (1411) . the upper bound for -F'Hoo 
from (]4"2l and the hypothesis ||fe||t»,i < P to the estimate (1221) with respect to the 
pairs (a', k') and (a, k). With = j( x ' ,v' )(t) , we obtain 

I A; — k\(y(t), y(t), w) du y (w) dt 



Sy(t)M 



< e(l + 2diam(M)pw„_ ie diam(M)s ) exp(2 diam(M) (e 
= eCi, say 



3 diam(M)S 



Co 



Finally, integrating the formula above in (x ,t> ) G T„ with the measure dfj,(x' ,v' ), 
we get 

Hifc-Jfe'lli < eVol(9M)w n _iCi. 
Theorem O holds now with C = max{Vol(9Af)w n _iCi, e diam ^ s /c }. 

7. Stability of the equivalence classes in two dimensions 

We prove here Theorem 14 .2\ making use of the results of [14"] . 

Let (a,k), (a,k) G V^, p be given with ||^4 — A\\* = e. As before, define the pair 
(a',k') in the equivalence class of (a, A;) by ( 13 7p and ( I38p . Then the corresponding 
albedo operator A' = A and, thus, 

(43) \\A' - A\\* = \\A- A\\* = e, and \\0 - p')\ sin^HU < e. 

The estimate (HO]) holds as in the case n > 3. Now 



(44) 

<^(x', f'; 

From (154"]) and (155]) . with y(s) = 7( x „)(s — r_(x, u)) 



--(x,v) 



(a - a)(y,y) ds 



< 



< 



~-(x,v) 



(\a! — a\ + \a — a|)(y,y) ds 

ds + £e diam(A./)S < 2£e diam(M) S . 



-{x,v) e diam(M)S 



T(y,y) 
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the same estimate holds for the second exponent in ( PHI) . Thus, from the definition 
QUE]) and flU we obtain 

(45) < exp{4ee diam ( M ) s } Hfc'IU < pexp{45e diam ( A/ ) E }. 
<p{x',v') 

Let 

Ex(y,w',w) := E(^ v>w > ) (-T-(y,w')),y)E(y,^ M (r + (y,w))) 

be the total attenuation along the broken geodesic due to one scattering at (y, w/, u>) G 
S 2 M. Then gl]) and g2J say 

(46) iw', iw)| > exp(-2(ee diam(M)s + diam(M)S)) = d, say, and 

(47) 1 1 - E[\\ < 2ee diam(M)s . 

The terms in the expansion of the albedo kernel in Proposition 13.21 are the 
traces of distributions defined on SM x T_; see [13]. The <f>j are the kernels of the 
operators J, KJ and (I — K)~ l K 2 J (j = 0, 1, 2, respectively) where 

Jf-(X,V) = E(^ x>v) (-T-(x,v)),x)f-(% jV )(T-(x,v))), 

fT— (x,v) 

Kf(x,v) = / S(x,7( a . iW )(t-r_(x,i;)))ri/(7( a , |1 ,)(t-r_(x,i;)))dt, with 
Jo 

T 1 f(x,v)= / k(x,v',v)f(x,v , )du x (i/). 

JS X M 

Let 7 be the trace operator on L°°(r + ), which is shown in [T3] to be well-defined. 

Let (sc, f , x', v') G r + x r_ be such that the geodesies 7( Xi „)(-) and J( x ',v')(-) intersect 
at (y,w',w) G S 2 M. By Proposition 13.21 above. 

E[(k - k') = (E[ - E x )k + (Eik - E[k') 

= (E[-E 1 )~k + CP- sin^l + (702 - 7<fe)| sin^|, 

and so by flU, (USD, flUD and (J37J), 

(48) Ci|jfe - k'\ < 25e diam(M)s p + e + | 7 0' 2 - 70 2 || sin# 
Now 

|70 2 - 70 2 | = l{I - KT'K'^'o ~ l(I ~ K)~ l K 2 (j)Q 

= 7(7 - K'Y 1 (K' 2 ^ - K 2 4> ) + 7(1 - K)-\K' - K)(I - K')~ y K 2 §q 
= 7(1 - K'Y 1 {K' 2 - k 2 )<f)' + 7(1 - K'r'K 2 ^ - O ) 

(49) + 7(7 - R)-\K' -K)(I - K')- x K 2 fo. 

Lemma 9 of [H] estimates the first of these terms: 

(50) || 7 (J - K')-\K' 2 - K 2 )^ < C 2 \\k' - fclUdlfc'IU + fk\\oo)(l - log | sinVI). 
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For the second we appeal to Proposition 7, and its proof, in [13]. Instead of using 
the estimate ||-E||oo < 1, we use (|42|) ; together with the fact that (J — K')' 1 preserves 
L°° , we readily obtain 

(51) || 7 (I - KT'K 2 ^ ~ to) < C 3 \\k\\U& -E\\oo< C 3 p 2 2ee d ^ M ^. 
The final term in ( l4"9j) is estimated in [1~4"| Lemma 10]: 

(52) || T (/ - Ry\K> -K)(I- KT'K'Moo < C A \\k' - fc|Up 2 (l + p). 
Combining ( 1481) . ( 150]) . (151~]) and ( 152]) . for a new constant C, we obtain 

W-Hoo^Ce + Cp\\k'-k\\ 00 
and so if p < 1/C, we obtain the final estimate 

C 
1 — Cp 

Theorem 14.21 is now proven for a constant which is the maximum of the C above 
and exp(diam(M)S)/c (see (HOI)). 
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